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We study the interplay between the chiral and the deconfinement transitions, both at high tem- 
perature and high quark chemical potential, by a non local Nambu-Jona Lasinio model with the 
Polyakov loop in the mean field approximation and requiring neutrality of the ground state. We 
consider three forms of the effective potential of the Polyakov loop: two of them with a fixed de- 
confinement scale, cases I and II, and the third one with a dependent scale, case III. In the cases 
I and II, at high chemical potential n and low temperature T the main contribution to the free 
energy is due to the Z(3)-neutral three-quark states, mimicking the quarkyonic phase of the large 
Nc phase diagram. On the other hand in the case III the quarkyonic window is shrunk to a small 
region. Finally we comment on the relations of these results to lattice studies and on possible com- 
mon prospects. We also briefly comment on the coexistence of quarkyonic and color superconductive 
phases. 

PACS numbers: 12.38.Aw, ll.lO.Wx, ll.30.Rd, 12.38.Gc 



I. INTRODUCTION 



Color confinement and chiral symmetry breaking are some of the most intriguing topics in modern theoretical 
physics. Quantum Chromodynamics (QCD) is believed to be the ultimate theory describing strong interactions. 
Nowadays it is accepted that the main ground state properties of QCD can be described in terms of non pcrturbative 
spontaneous breaking and/or restoring of some of the global symmetries of the QCD lagrangian. 

Unfortunately solving QCD in its non perturbative regime is a hard task. At zero and small quark chemical potential 
fi lattice calculations are a good tool to derive the equation of state of QCD matter, the transition temperatures and so 
on starting from the first principles, see for example and references therein. Several approximation methods 

are available to overcome the sign problem of the fermion determinant with three colors at finite fj. (see Refs. 0, H 0| 
for reviews on the sign problem) : small- expansion [1, [13] , reweighting tecniques [l^ , density of the states 
methods [T^ and analytic continuation to imaginary chemical potential [fl [isl. [la. [TtI. Tlg| . 

Besides lattice calculations there exist effective descriptions of QCD. Among them Nambu-Jona Lasinio (NJL in 
the following) models [I0| are very popular, see [23| for reviews. They are based on the observation that several 
properties of the QCD ground state are related to the spontaneous breaking of some of the global symmetries of the 
QCD lagrangian. Therefore one hopes that by a model that has the same global symmetry breaking of QCD one can 
capture the essential physics of QCD itself. 

In recent years it has been argued that the NJL model, which does not contain gluons, can be improved by adding 
a non linear term to the lagrangian which describes the dynamics of the traced Polyakov loop [2l| , and an interaction 
term of the Polyakov loop with the quarks. The resultin g rn odel is called the PNJL model, introduced in Refs. (2^.[23| 
and extensively studied in [H [H [ll, [13, [H, [H, M, [MIH [H, [H, M, M, [13, M, M, \M, \M, [13, [H, [3 ■ In the PNJL 
model one assumes that a homogeneous euclidean temporal background gluon field couples to the quarks via the 
covariant derivative of QCD. This coupling gives rise to the interplay between the chiral condensate and the Polyakov 
loop. Even if it is very simple, the PNJL model turned out to be a powerful tool which allows to compute several 
quantities that can be computed on the lattice as well. The agreement with existing lattice data is satisfactory p3 |. 

One of the exciting characteristics of the PNJL model is the statistical confinement of quarks at low temperature ^ . 
In a few words this means that at small temperature and small chemical potential the contribution to the free 
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energy, fl, of the states with one and two quarks are suppressed, and the leading contribution to fl arises from the 
thermal excitations of colorless three quark states. This property is related to the small value of the expectation value 
of Polyakov loop which is found in the self-consistent calculations within the PNJL model in the aforementioned 
conditions of temperature and chemical potential. It has been recently argued by Fukushima that the statistical 
confinement property of the PNJL model persists even at high chemical potential [ 33i] . This result is in agreement 
with the phase diagram of QCD obtained in the large number of colors Nc approximation 0,14^, see also Refs. [47Ll48l| 
for recent related studies. Inspired by Ref. ^4511 Fukushima has suggested to interpret the statistical confined phase of 
the PNJL model at high quark chemical potential as the quarkyonic state found in [4^ . 

In this work we investigate on the ground state of the electrically neutral two flavor PNJL model, focusing on 
its possibl e q uarkyonic structure at high /j, and low T. We use a non local four fermion interaction instead of the 
local one [l9|. The local NJL model is usually regularized by means of an ultraviolet sharp cutoff, which amounts 
to artificially cutoff the quark momenta that are larger than the cutoff itself. Thus the extensions of the model to 
temperatures and/or chemical potentials of the order of the cutoff are quite dubious. However if one introduces a 
non local interaction, which corresponds to the multiplication of the NJL coupling by a momentum dependent form 
factor f{p), and requires that the form factor satisfies the asymptotic freedom property of QCD f{p oo) — 0, then 
all of the momentum integrals are convergent and the model is consistent at any value of temperature and chemical 
potential. In this paper we use one specific form of the form factor. Although the choice of a different functional form 
for f{p) can lead to different quantitative results (mainly the shift of the critical points) we believe that our picture 
should not be modified qualitatively. We consider the logarithmic form of the Polyakov loop effective potential U 
suggested by Ratti, Roessner and Weise in Ref. [2^; moreover we investigate on the effects of a dependence oiU on 
the quark chemical potential as well as on the number of flavors as suggested in Ref. [l^. We compare the phase 
diagrams obtained in the cases in which we do not consider (cases I and II) and do consider (case III) the ^-dependence 
of U. Cases I and II differ for the value of the deconfinement scale in the Polyakov loop effective potential. 

We find that the phase diagram in the two cases (I and II on one side. III on the other side) differ even qualitatively. In 
particular, on one hand in the cases I and II we confirm the results of Fukushima [s^l and strengthen his interpretation 
of the high chemical potential/small temperature state of the PNJL model as the quarkyonic matter of the large Nc 
phase diagram. On the other hand in case III we find that the quarkyonic-like window found in the cases I and II is 
shrunk and becomes a small region in the fJ. — T plane in the case III, opening a wide room for the dcconfined quark 
matter of the pure NJL model. 

The plan of the paper is as follows. In Section II we sketch the formalism. In Section III we discuss our results. 
Finally in Section IV we draw our conclusions. 

II. THERMODYNAMIC POTENTIAL WITH A NON LOCAL FOUR FERMION INTERACTION 

The Lagrangian density of the two flavor PNJL model is given by [H, |4l| 

£' = e(z7^a^)e + {ijf.D^' + ^^70 - m) + A - T] . (1) 

In the above equation e denotes the electron fleld; ip is the quark spinor with Dirac, color and flavor indices (implicitly 
summed), m corresponds to the bare quark mass matrix; we assume from the very beginning — rrid- The covariant 
derivative is deflned as usual as Z?^ = d^~iAf^. The gluon background field = 6ofiAo is supposed to be homogeneous 
and static, with Aq = gAgTa and Ta, a — 1, ... ,8 being the SU{3) color generators with the normalization condition 
Tr[Ta,T;,] = 6ab- Finally /i is the chemical mean quark chemical potential, related to the conserved baryon number. 

In Eq. ([T]) <&, $ correspond to the normalized traced Polyakov loop L and its hermitian conjugate respectively, 
$ = TtW/Nc, # = TiW''/Nc, with 



W exp ^ J A^dr^ = exp (ipA^) , A^ = iAo 



(2) 



and /? = 1/T. $ is a color singlet but it has a Z{3) charge [21j, where Z{3) is the center of the color group SU{3); 
thus if $ 7^ in the ground state then the Z{3) symmetry is spontaneously broken. The term U[^,^,T] is the 
effective potential for the traced Polyakov loop; in absence of dynamical quarks it is built in order to reproduce the 
pure glue lattice data of QCD, namely thermodynamical quantities (pressure, entropy and energy density) and the 
deconfinement temperature of heavy (non-dynamical) quarks, T = 270 MeV. Several forms of this potential have 
been suggested in the literature, see for example In this paper we adopt the following logarithmic 

form I25l. 



^2^$$ + b{T) log [1 - 6$$ + 4($3 + $^) - 3($$)^] 



(3) 
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with 



Numerical values of the coefficients are as follows [25|: 

ao = 3.51 , ai = -2.47 , as = 15.2 , 63 = -1.75 . (5) 

If dynamical quarks were not present then one should chose Tq = 270 MeV in order to reproduce the deconfinement 
transition at T = 270 of the pure gauge theory [H, [U, [2^. In presence of quarks Tq might get a dependence on the 
number of active flavors as well as on the quark chemical potential [23,[2^. Inspired by Refs. [23l [2^. in this paper 
we consider three cases: 

To = 208 MeV , Case I , (6) 
To = 270 MeV , Case II , (7) 
foi^i) = T^e-i/"«'=('^^ , Case III . (8) 

Case II corresponds to the deconfinement temperature in the pure glue theory; the parameters in the cases I and III 
have been evaluated in Ref. [2^ on the basis of hard dense and hard thermal loop approximations to QCD. In the 
equation corresponding to Case III we have set 

ao = 0.304 , Tr = 1770 MeV , (9) 

and 

with N j = 2 and A^c = 3. At ^ = we have To(/^ = 0) = 208 MeV as case I; for comparison, at = 500 MeV the 
deconfinement scale is given by To(/x = 0) = 19 MeV. 

In Eq. ([1]) £4 represents the lagrangian density for the four fermion interaction. If we define S4 — J d^xC^ as the 
interaction action then in the local version of the NJL model one has 



S^^GJ d^x [(^V)' + ii^^ToTi^f] ■ (11) 
In the non local version of the NJL model the contact term Eq. pT|) is replaced by [13, 113, HO, [U, [H, [H, [13] 

S^^gJ d^x [{q{x)q{x))' + iq{x)i^5rq{x))'] , (12) 
where the dressed quark field is defined as 

qix) = / d^y Fix - y)i^{y) , (13) 



and F{r) is a form factor whose Fourier transform f{p) satisfies the constraint f{p) — > for p — > 00, p being the 
3-momentum. In this paper we follow Ref. [30j and use the Lorentzian form factor, 

f{p) = , ^ ^^^^^ . (14) 

v/1 + (J3/A)2" 

In the above equation A = 684.2 MeV and a = 10. Moreover we use m = 4.46 MeV and G = 2.33/ A'^ By these 
numerical values we reproduce the pion decay constant /^r = 92.3 MeV and the pion mass m-^ = 135 MeV, as well 
as the chiral condensate (uu) = —(256.2 MeV)'^. Although the choice of a different form factor will lead to different 
critical temperatures and/or chemical potentials, it is quite reasonable that the qualitative picture that we draw in 
this work is insensitive to the specific form of f{p)- 

As explained in the Introduction we are interested to the ground state of the model specified by the Lagrangian in 
Eq. Il]), at each value of the temperature T and the chemical potential /i, corresponding to a vanishing total electric 
charge. In order to build the neutral ground state we use the standard grand canonical ensemble formalism, adding 
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to Eq. JT]) the term ^jlqNq, fj,Q being the chemical potential (i.e. Lagrange multiplier) for the total charge Nq, and 
requiring stationarity of the thermodynamic potential with respect to variations of /xg, which is equivalent to the 
requirement < A^g in the ground state. This amounts to write the lagrangian C in the gran canonical ensemble 
C = C' + ^1QNQ as i4l| 

C = e(n^9^ + ^lelo)e + i^ [il^,D^' +fi-io-m)^ + G [{iji^f + {i^i-f^f^f'\ - U[<^, $, T] , (15) 
where = ~I^Q and the quark chemical potential matrix fi is defined in flavor-color space as 

where Ic denotes identity matrix in color space. At /ig 7^ a difference of chemical potential between up and down 
quarks, 5^, — /J.2/2, arises. 

In this paper we work in the mean field approximation. Because of 7^ a pion condensation might occur 
in the ground state [ssl ]. In order to study simultaneously chiral symmetry breaking and pion condensation we 
assume that in the ground state the expectation values, real and independent on x, for the following operators may 
develop [H, [M [H, [sl, [13 , 

cr G{q{x)q{x)) , tt = G(g(a;)«75Ti9(a:)) . (17) 

In the above equation a summation over flavor and color is understood. We have assumed that the pion condensate 
aligns along the ri direction in flavor space. This choice is not restrictive. As a matter of fact we should allow for 
independent condensation both in tt"*" and in 7r~ channels [3l|: 

^± = G(^z75r±^) = ^e±'^ , (18) 

with T± = (ti ±r2)/V2; but the thermodynamical potential is not dependent on the phase 0, therefore we can assume 
9 = Q which leaves us with 7r+ = tt^ = 7r/V2 and introduce only one condensate, specified in Eq. (|17p. 

In what follows we consider the system at finite temperature T in the volume V . This implies that the space-time 
integral is / d'^x — dr J d^x with /3 — 1/T. In the mean field approximation the PNJL action reads 



S ^ J d^x [e{i-ff,^^' + Aie7o)e + {h^.D^' + fljo) V'] 
+2(7 / d'^x q{x)q{x) + 27r / d'^x q{x)ij5Tiq{x) 



^2 



-PV—^ PVU[^,<f,T], (19) 

where V is the quantization volume and (3 = l/T. In momentum space one has 

^ ^ J (2^ [^(^mP^ + Me7o)e + (7mP^ " 7/^^^ - A7o) V'] 
[2o- ^p{p)^p{p) + 27r V(p)i75nV'(p)] 

- PV^—^ - $, T\ , (20) 

with yip ~ gAljTa- We introduce the mean field momentum dependent constituent quark mass M{p) and renormalized 
pion condensate N{p): 

M{p) = m - 2<jf{p) , N = -2Tif{p) . (21) 

The thermodynamical potential per unit volume in the mean field approximation can be obtained by integration 
over the fermion fields in the partition function of the model, see for example Ref. (56l |. 

/ u4 ^2j.2 7^2j.4X _ ^2 ^^2 



127r2 6 180 / L ' ' J Q 



^E/I^THog^:^, ,22, 



5 



where the sum is over fermion Matsubara frequencies aj„ = 7rT(2n + 1), and the trace is over Dirac, flavor and color 
indices. The inverse quark propagator is defined as 



S ^{iuj„,p) = 



(23) 



Performing the trace and the sum over Matsubara frequencies we have the effective potential for $, a and tt, namely 



n = - 



127r2 



-2T 
-2T 

-2T 



(2^ 
d^p 
(2^ 
_fp_ 
(2^ 
d^p 
(2^ 
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log 



180 



-Z^[$,I>,T] 



G 



f2^ 



[£;+ + s_ - 2p] 



log 
log 
log 



where 



7V2 



(24) 
(25) 



and Ep — -y/p^ + A'P{p). In Eq. (|24|) the integral of 2p is an irrelevant constant that we subtract in order to make the 
thermodynamical potential finite at each value of temperature and chemical potential. The ground state of the model 
is defined by the values of cr, tt, <&, $ that minimize f2 and that have a vanishing total charge; the latter condition is 
equivalent to the requirement 







In this paper we use the convenient Polyakov gauge, 

1, 



$ = -Tr 
3 



(26) 



(27) 



with <^3, (/jg real parameters. It has been widely discussed in Ref. (25| that in the mean field approximation and with 
the choice of the effective potential U given by Eq. ^ one has (<&) = ($) for any value of T and /i, and the solution 
($) ^ ($) at finite /i is due to quantum fluctuations. Since in this paper we consider only the mean field approximation 
we chose $ = $ in the calculations. This choice implies (/)8 = thus we are left with only one parameter (j)^ = 0. 

Before closing this section we write the dispersion laws of the quasi-particles in the Polyakov gauge, defined as the 
poles of the quark propagator given by Eq. P5|) : 

(28) 
(29) 
(30) 

In the previous equations m, d correspond to up and down quarks, r, g and b to the colors red, green and blue; the 
upper (lower) sign multiplying ^ and <j) correspond to quarks (antiquarks) . 





= =P/i ±i(f) - 


^E+ , 


Edr 


= ± - 


f 


E 




^E+ , 


Edg 








Eub — ~ 




Edb 


= TM + ^- 





III. SUSCEPTIBILITIES IN THE PNJL MODEL 



In order to study the landscape of the phases of the PNJL model we introduce the susceptibility matrix. Suscepti- 
bilities are useful to identify phase transitions since they are proportional to the fluctuations of the order parameters 
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around their mean field values, which usually are enhanced near a phase transition. We follow closely Ref. ^3d\ for the 
formalism settings. The first step is the definition of the dimensionless curvature matrix of the free energy around its 
global minima C [H, [s^l , 

(Cmm Cm'S> \ 
Cm* C$$ C$|, . (31) 
Cm* C$(j, Cq,q, J 

In the above equation the diagonal entries are defined as 

Cmm - XaM2 ' (32) 

** ~ A3 ' ** ~ A3 ' y-^-^' 

with f3 = 1/T and A is the mass scale defining the form factor Eq. (fT4|) . fl is defined in Eq. ((24|) . In what follows we 
denote by M the constituent quark mass computed at p = 0, which is a function of fi and T. The off diagonal entries 
are given by 

_ P d'n _ (3 d'n 

" A3 9ci>9$ ' ^-^^^ 

the derivatives are computed at the global minimum of Q. Notice that the proper definition of the curvature matrix 
requires that we put $ = namely the mean field solution, only after differentiation. 

The susceptibility matrix x is computed as the inverse of the curvature matrix C. We have 

(Xmm Xm* Xa/* 
XM* X** X** I ■ (36) 
Xm* X** X** 

Here xmm, X** x*$ denote respectively the dimensionless susceptibilities of the constituent quark mass, of the 
Polyakov loop and of its complex conjugate. We also introduce the average susceptibility 

X = ^ (X** + X** + 2x**) • (37) 



IV. RESULTS AND DISCUSSION 



In this Section we sketch our results. Firstly we discuss the set of parameters corresponding to the case I which 
corresponds to Tq — 208 MeV. Case II is qualitatively similar to case I, therefore after the discussion of the results 
obtained in the latter case we briefly show the results corresponding to the former case. Finally we compare both 
qualitatively and quantitatively the cases I and III. We find that the phase structures of the models corresponding to 
cases I and III are quite different. 



A. Case I: masses, Polyakov loop and quarkyonic matter 

In the upper panel of Fig. [T] we plot the constituent quark mass at p = 0, the expectation value of the traced 
Polyakov and the electron chemical potential as a function of the temperature, computed at /i = (left) and = 300 
MeV (right). AIq denotes the constituent quark mass at p = 0, /Li = 0, /^e = and T — 0, Mq = 335 MeV. The pion 
condensate N is not shown since we find = once electrical neutrality has been imposed. The latter result is in 
agreement with what we have found in our previous work, see Ref. [4l'|, where we have considered the local version of 
the neutral two flavor PNJL model. Even if we have shown results only for two values of the quark chemical potential, 
we have explicitly verified that N vanishes in the whole range of chemical potentials and temperatures considered in 
this work, namely < < 500 MeV and < T < 250 MeV. 
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FIG. 1: Upper panel: constituent quark mass at p = 0, pion condensate at p = and Polyakov loop as a function of temperature, 
computed at ^ — (left panel) and /j. = 353.5 MeV (right panel). Mq denotes the constituent quark mass at p = , /x = 0, 
fj,e = and T — 0, Mo = 335 MeV. Lower panel left: susceptibilities at ^ = as a function of temperature. Lower panel right: 
susceptibilities at /i = 300 MeV as a function of temperature. Solid line: xmm- Dashed line: Xsi.- Dot-dashed line: x- 



The expectation value of the Polyakov loop at /i = is consistent with zero up to temperatures of the order of 100 
MeV. ^ It raises as the temperature is increased becoming of the order of 1 for temperatures close to 250 MeV. This 
behavior signals a crossover from a low temperature phase with an unbroken Z(3) symmetry, to a high temperature 
phase with Z{3) symmetry spontaneously broken. The behavior of as a function of the temperature is observed even 
at higher values of /i, see for example the right upper panel of Fig. [1] We call such a crossover as the Z{3) crossover 
throughout this paper. 

In the lower panel of Fig. [1] we plot three of the susceptibilities defined in the previous Section, namely xmm (solid 
line), x$<E) (dashed line) and x (dot-dashed line), as a function of temperature at = (left) and fi = 300 MeV (right). 
In this work we identify the chiral crossover temperature with the temperature where xm m is maximum. In the same 
way and following Ref. [30| we define the Z{3) crossover temperature as the one corresponding to the maximum of x- 

We wish to investigate on the spontaneous breaking of the Z{3) symmetry in the neutral PNJL model as the quark 
chemical potential is increased at a fixed low temperature. To this end we plot in Fig. [2] the constituent quark mass 
at J) = 0, the expectation value of the traced Polyakov and the electron chemical potential as a function of the quark 
chemical potential n, computed at T = 20 MeV (left panel). Mq denotes the constituent quark mass a,t p = 0, fi = 
and T = 0, Mq = 335 MeV. Again we do not show the pion condensate since it turns out to vanish in the neutral 
phase. At low temperatures we find a first order chiral transition at /i ~ 353 MeV, in agreement with our previous 
analysis [4l|. In correspondence of the chiral restoration the expectation value of the Polyakov loop has a sudden 
jump. Nevertheless its value remains much smaller than one even if fi is increased to 500 MeV, where $ w 0.04. For 
comparison we show the same quantities at T = 130 MeV in the right panel. 

We now focus on the low temperature regime, therefore we refer to the left panel of Fig. [2l In this case we can 
not identify the jump of $ as the Z{3) crossover. Instead the discontinuity of <!> is simply due to the coupling of the 



^ <J> can not be exactly zero because dynamical quarks break the Z{3) symmetry explicitly; nevertheless <J> turns out to be very small, 
signaling that the center symmetry is broken only softly. 
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FIG. 2: Upper panel: constituent quark mass at p = 0, Polyakov loop and electron chemical potential as a function of the quark 
chemical potential ^, computed at T = 20 MeV (left) and T = 130 MeV (right). Mo denotes the constituent quark mass at 
p = and T = 0, Mo = 335 MeV. In both cases A'^ = and it is not shown. Lower panel: susceptibilities at T = 20 MeV (left) 
and T = 130 MeV (right) as a function of the quark chemical potential. Solid line: Xmm- Dashed line: Dot-dashed line: X- 



Polyakov loop with the chiral condensate. This is confirmed by the calculation of the Polyakov loop susceptibilities, 
see the lower panel of Fig. [2l At T = 20 MeV, in correspondence of the jump of the constituent quark mass, the chiral 
susceptibility has a pronounced peak. On the other hand the Polyakov loop susceptibilities are very smooth functions 
of ^ with a small cusp in correspondence of the chiral transition, signaling the absence of a phase transition (as well 
as of a crossover). For comparison we show the same quantities at T = 130 MeV in the right panel. 

Our results can be interpreted by assuming that at low temperatures the Z[i) symmetry is not spontaneously 
broken, both at low and at high chemical potentials. The non zero value of $ can be related to the existence of 
dynamical quarks in the system, that break explicitly the center symmetry. The fact that $ 1 means that in 
the ground state colored quarks are suppressed (they have a finite Z(3)-charge), and the main contribution to the 
free energy is due to the .Z(3)-invariant multi-quark states, that are states with a zero Z'(3)-charge. This point can 
be clarified by studying the thermal population of the quasi-quarks excitations at low temperature. To this end we 
compute the quark number density rig, 

as a function of the chemical potential at fixed temperature. The result is shown in Fig. [3] Evaluation of the derivative 
of Vl defined in Eq. ([M]) leads to the expression 



I p dp 



9+_ 9^ 
/+- /— 



g+- 



g-+ 



/++ /- 



where we have introduced the functions 



/±± = 1 + Sc&e-'^^^ii^) + 3<i>e-2/3(^=±±A') 



-3/3(_B±±m) 



3/3(i5±±p) 



(39) 



(40) 
(41) 
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FIG. 3: Left panel: quark number densities, in units of the nuclear saturation density po — 0.16 fm~^, as a function of the 
quark chemical potential at T = 20 MeV (dashed line) and T = 130 MeV (dot-dashed line). Right panel: dimensionless 
quark number susceptibilities as a function of the quark chemical potential p at T = 2Q MeV (dashed line) and T — 130 MeV 
(dot-dashed line). 



and E± are defined in Eq. ([25]). The addenda in the r.h.s. of Eq. ((39|) correspond respectively to up quarks, down 
quarks, up antiquarks and down antiquarks. If we put by hand $ = 1 in Eq. (|39|) we recover the usual expression of 
the NJL model. 



3 f°° 
nq,NJL = — / P^dp 
Jq 



1 



1 



1 



1 



(42) 



where the 3 overall counts the number of colors. Eq. (|42p is the number density of a free fermion gas; it shows that in 
the zero temperature limit and fi > M, M denoting the constituent quark mass, the ground state of the NJL model is 
made of Fermi spheres of red, green and blue quarks. Moreover at small but non vanishing temperatures the thermal 
excitations over the Fermi spheres are still quarks. 

Now we compare Eq. (j42|) with the analogous result of the PNJL model. At low temperature we have <& <C 1 
therefore for a rough analysis we can put $ = in Eq. ([55]) . We are left with the expression: 



T-q.PNJL — 



p^dp 



1 



1 



1 



1 



(43) 



The above equation is valid for every value of fj,. In the limit T ^ and for n > M, with M the constituent quark 
mass, it gives the equation obtained in the NJL model, that is a ground state of Fermi spheres of red, green and blue 
quarks at the chemical potential /i. If we introduce a small temperature then the thermal excitations are not quarks 
but the Z(3) symmetric three quark states, that is states made of one red quark, one green quark and one blue quark. 
This is clear from the above Eq. by looking at the arguments of the exponentials in the four addenda. Each of 
the addenda corresponds to the occupation number of fermions with energy given by 3E± — 3/i which is exactly the 
energy of the lightest Z{3) symmetric state, namely (see Eqs. (P5)) - (P(I)) ) 



Er, 



■ Eu-L 



3E± - 3n 



(44) 



the sign depending on the flavor we consider (_E+ corresponds to up quarks, i?_ to down quarks). The same result 
holds for antiquarks, simply by replacing fi — > — /i. The combination (j44p is exactly the argument of the exponentials 
in Eq. dSl). 

To summarize: for the parametrization I the ground state of PNJL quark matter in the regime of low temperature 
T <^ M and fi > M is made of Fermi spheres of quarks, and the thermal excitations above the aforementioned Fermi 
spheres are the three quark states, neutral with respect to Z{S). 

For completeness, on the right panel of Fig. Owe plot the dimensionless quark number susceptibilities, Xg, defined 

as 



Xq 



1 d^n 



(45) 



where A is the form factor momentum scale in Eq. (HH), and il is the PNJL free energy given by Eq. (I24p . 
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FIG. 4: Phase diagram of the neutral two flavor PNJL model. The dot-dashed line corresponds to the chiral crossover; the bold 
solid line is the flrst order transition. The thin solid line denotes the deconfinement crossover. 



B. Case I: phase diagram in the /i T plane 

In Fig. [5] we summarize the phase diagram of the model in the fJ, — T plane with the parametrization I. The thin 
line corresponds to the chiral crossover; the thick line is the first order chiral transition. We identify the peaks (or the 
local maxima) in the susceptibilities with the phase transitions. In particular, the chiral crossover is related to the 
peak of XMM in Eq. ([36| ; on the other hand following Ref. [s^] we identify the peak of the average susceptibility x 
defined in Eq. ((37)) with the Polyakov loop crossover. 

From the qualitative point of view the phase diagram does not differ from our previous result [4ll | obtained in the 
sharp cutoff regularization scheme. The chiral crossover at /i = is located at Tc — 215 MeV, to be compared with 
our previous work [4l| Tc = 206 MeV. The critical end point is only slightly shifted: in this work we find 

{fiE,TE)~ {350,55) MeV ; (46) 

this result has to be compared with [4l| 

{Pe,Te) w (340, 80) MeV , sharp cutoff. (47) 

Finally at T = wc find that the chiral crossover occurs at /i = 370 MeV, while in our previous work we have found 
H = 350 MeV. 

We now discuss the Polyakov loop crossover line, corresponding to the thin solid line in Fig.[4l At small values of the 
quark chemical potential the peaks of the averaged susceptibility are well pronounced, see for example Fig. [TJ As /i is 
increased the peaks of x s-s well as of the diagonal , x^^ ^-iid off-diagonal x^<^ susceptibilities are broadened and 
the crossover is dilute over a wide interval of temperatures, see the right panel in Fig. [T] In the window of chemical 
potential studied in this paper, < /i < 500 MeV, we are still able to observe maxima of x (as well as for the other 
susceptibilities) as a function of the temperature at a fixed value of /z; the width of the maxima increases as /i is 
increased. Therefore we expect that at high values of /i and T the peaks of x will be very dilute, meaning that the 
crossover disappears in the model under consideration. This result changes if we consider /i dependent coefficients of 
the Polyakov loop effective potential as we discuss later. 

We finally notice that our results for the Z{3) crossover is in qualitative agreement with the results obtained in 
Ref. [soj, where the authors study the phase diagram and the susceptibilities of the PNJL model with quarks at the 
same chemical potential, and with a polynomial form of the Polyakov loop effective potential U. This suggests that 
the Z{S) crossover is not mainly governed by the specific form of U or by electrical neutrality, but by the assumption 
that the deconfinement scale Tq in U is kept independent on fi in this calculation. 
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FIG. 5: Upper panel: constituent quark mass at p = 0, Polyakov loop and electron chemical potential as a function of the quark 
chemical potential ^, computed at T = 20 MeV in the case III (left) and case I (right, shown for comparison with case III; it 
is the same plot shown in Fig. [2}. Mo denotes the constituent quark mass at p = 0, = and T = 0, Mo — 335 MeV. In both 
cases = and it is not shown. Lower panel: susceptibilities at T = 20 MeV in the case III (left) and case II (right). Solid 
line: Xmm- Dashed line: Dot-dashed line: x- 



C. Case II: critical points 

From the qualitative point of view the case with Tq — 270 MeV does not differ from the previously analyzed case 
II. Therefore we simply give the coordinates of the critical points obtained in this case. At fj, — we find the chiral 
crossover at T = 219 MeV and the Z{3) crossover at T = 211 MeV. The critical end point coordinates are 

ifiE, Te) « (336, 103) MeV , Tq = 270 MeV . (48) 



D. Case III: critical points and phase structure 



We now discuss the results obtained in the case III in which we assume both a /x and a Nf dependence of the 
parameter Tq of the Polyakov loop potential, see Eq. ([8]). Our main goal is to emphasize the differences between case 
III and case I. The main difference arises at low temperature and high chemical potential, so we focus on this regime. 
In Fig. [5] we plot in the left panel the constituent quark mass at p = and the expectation value of $ as a function 
of /i at T = 20 MeV, with the related susceptibilities, for the case III, and compare these results with those obtained 
in the case I at the same temperature (right panel). We have verified that qualitatively the picture does not change 
if we lower the temperature to the order of one MeV. 

At /i = the critical temperatures are equal to those computed in case I (simply because To(/^ = 0) = 208 MeV). 
Moreover the coordinates of the critical end point are 

{^iE, Te) = (339, 53) MeV , To = %{^x) . (49) 

The data on $ corresponding to the parametrization III show that the case Tq — TQ{fi) is quite different from the 
case To = 208 MeV. In the case III (left panel) in correspondence of the chiral transition at /i = /ic ~ 350 MeV the 
Polyakov loop has a net jump from <E> <C 1 at /x = /^c — 0^ to a definitely non zero value $ w 0.3 at fi — fic + Q'^ . Since 
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the contribution of the one and two quark states {Z{3) charges) to the free energy is multiphed by 3$, see Eq. (|24p . 
and in the present case 3$ is of the order of unity, the weight of the Z{3) charges in the free energy is the same of the 
weight of the three quark states. This behavior is different from what we have found in the case of Tq = 208 MeV. 
The similarity between the two cases is partially recovered if we consider temperatures of the order of one MeV; in 
this case we find a narrow window in fi where 3$ is of the order of 0.1, revealing a ground state in which the leading 
contribution to the free energy comes from the thermal excitations of Z{3) neutral states. We discuss this point in 
more detail in the following Section. Finally the analysis of the peaks of the susceptibilities Xi<^ ^nd x (lower left 
panel in Fig. [5]) reveals that the Z{3) crossover occurs at /Lt « 460 MeV. 

V. COMPARISON BETWEEN THE TWO SCENARIOS 

In this Section we compare the qualitative picture that arises from the study of the phase diagram of the neutral 
PNJL model within two scenarios: the first one corresponds to keeping an independent Tq, case I; the second one 
corresponds to keeping a /^-dependent Tq, case III. 

The results that we have discussed in the previous Sections show that the phase diagram of the PNJL model in 
the case I at low temperatures is similar to the phase diagram obtained in the large Nc approximation of QCD, see 
Refs. [i^ H^. At low temperatures the latter phase diagram consists of two regions: the first one at low values of fi, 
defined as the confined phase and characterized by $ = and a vanishing baryon density; the second one at large 
values of /i called quarkonia in which $ = but the baryon density is not vanishing. Finally at high temperature one 
finds the deconfined phase with $ 7^ and a non vanishing baryon density. In the quarkyonic phase the free energy 
is that of free quarks, but the thermal excitations are those of baryons. Our previous discussion and Figs. [2] and [3] 
show that this happens even in the PNJL model in the low temperature regime. Therefore the PNJL model with 
parametrization I approximately reproduces the large Nc phase diagram at low temperatures, if one interprets the 
state with $ ^ 1 at high /i with the quarkyonic phase of large Nc- This fact has been already noticed in a study 
of the three flavor model by Fukushima [33] where the author has suggested to identify the low temperature-high 
density ground state of the model as the quarkyonic phase of large Nc QCD. Our results strengthen this idea and 
thus suggest that the quarkyonic-like ground state of low temperature-high density PNJL model is not a peculiarity 
of the three flavor case, but it seems to be a characteristic of the PNJL model itself, as far as we do not include 
an explicit fj, dependence into the coefficients oi U (we discuss this case in a next Subsection). The main difference 
between large Nc and PNJL is that in the latter model one can excite one and two quark states (that is Z{3) charges) 
if the temperature is high enough. As a consequence, the deconfinement transition observed in the large Nc model at 
high temperature and high chemical potential is replaced in the present model by a smooth Z{3) crossover. 

In Fig. [5] we show a cartoon phase diagram of the neutral two flavor PNJL model and a comparison with that 
obtained in the large Nc approximation [4^ • The bold line denotes the chiral crossover as well as the chiral first order 
transition. The thin line corresponds to the deconfinement crossover. Both of these lines are the same which we have 
shown in Fig. 21 Since this is simply a cartoon we do not distinguish between the crossover (small ^) and first order 
transition (higher values of n). In the PNJL model the quark density does not vanish at any finite temperature, even 
if for small chemical potential Ug is very small at low temperature (see Fig. [3]). To compare the phase diagram of the 
PNJL model with that of the large Nc approximation we need a criterion to say if Uq is zero or not. Analogously to 
Rcf. 33] we identify the Uq crossover with the value of ^ corresponding to the inflection point of the quark density. We 
find that the Uq crossover defined in this way coincides with the chiral crossover as in Ref. [3^. Therefore the chiral 
crossover line in Fig. [5] represents the density crossover as well. In the chiral broken phase and at low temperature 
Uq ~ 0. On the other hand ^ in correspondence to the chiral symmetric phase at low temperature. At high 
temperature ^ both in the chiral broken and in the chiral restored phases. 

At low temperature we have $ « both on the left and on the right of the dashed line, see Figs. [T] and [H At 
low temperature the region with broken chiral symmetry has the same characteristics of the hadron phase found in 
Ref. [4^ ; on the other hand at low t emp erature the region on the right of the dashed line has the same characteristics of 
the quarkyonic phase found in Ref. [45|. For these reasons we have called the two regions hadronic-like and quarkyonic- 
like respectively. We stress that this analogy holds strictly speaking only at low temperature (for temperatures of the 
order of one hundred MeV ^ even in the chiral broken phase, see Fig. [3]). Finally at high temperature (above 
the Z{3) transition line) we have both $ of order of unity and Uq ^ 0. In analogy to the terminology of Ref. [4^ we 
call this region of the phase diagram the deconfined-like phase. 

We briefiy compare the results discussed above in relation with the case I with those obtained in the large Nc 
approximation and at T = in Ref. [48j . where the author discusses a gap in the spectrum of quarkyonic matter 
within a model. Such a gap is given by the pion mass, M^, which becomes larger as fi is increased. Even if the values of 
as a function of /j, computed in Ref. [i^ might differ from the non local PNJL ones, the calculations of carried 
out in Refs. [13, |4l| using the local NJL model show that the qualitative behavior of as a function of ^ is the 
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FIG. 6: Cartoon phase diagram of the neutral two flavor PNJL model and comparison with that obtained in the large Nc 
approximation 45]. The bold line denotes the chiral crossover as well as the chiral first order transition. The thin line corresponds 
to the deconfinement crossover. Both of these lines are the same which we have shown in Fig. [l] The n, crossover coincides with 
the chiral one. At low temperature we have "1> « and ~ in the chiral broken phase, in agreement with the hadronic phase 
of Ref. 14511 . At low temperature and in the chiral symmetric phase we find $ « and 7^ in agreement with the quarkyonic 
phase [45|. For these reasons we have called the two low temperature regions hadronic-like and quarkyomc-hke respectively. 




FIG. 7: Low temperature phase diagram of the PNJL model with parametrization III. Bold line denotes the chiral transition. 
Thin line corresponds to the Z(3) transition. The shaded region denoted by Q-L corresponds to the zone of the quarkyonic-like 
state of matter. 



same in the two models. Thus in the PNJL model we expect a large pion mass at large fi as well. However this mass 
does not correspond to the gap in the excitations spectrum in our model. As a matter of fact in the quarkyonic-like 
region of the phase diagram in Fig. [5] the three quark states can be excited; each quark has a constituent mass M{p) 
given by Eq. (|2ip and plotted in Figs. [l]and[2l hence the three quark state has a mass 3M{p) which at small quark 
momenta and large n is of the order of 10 MeV. Therefore in our case a gap in the spectrum still exists but it is given 
by the three quark state mass which is much lighter than M^. 

We now turn on the parametrization III. As discussed in the previous Section, the small chemical potential region 
of the phase diagram in case III is qualitatively similar to that obtained in the case I, therefore we focus on the low 
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temperature/large chemical potential region from now on. In Fig. [7] we draw the low temperature phase diagram 
of the PNJL model with parametrization III. The bold line denotes the chiral transition; the thin line corresponds 
to the Z{3) transition. As in the previous section the transition lines are computed by looking at the peaks of the 
chiral and x susceptibilities. The diagram in Fig. [7| should be compared with the analogous diagram obtained for the 
parametrization I which is shown in Fig. |4l The main effect of choosing the parameter Tq as a /i dependent one in the 
Polyakov loop potential is the lowering of the Z{3) transition line. Moreover the wide quarkyonic-like window in Fig. [6] 
is shrunk to a small region in Fig. [T] At low temperature it is enough to reach a chemical potential of the order of 
500 MeV to have $ « 1 and a net quark density; both these characteristics define the deconfined phase of Fig. [5] 4j|. 
Even if we have used the particular form Tb(/x) suggested in Ref. [1^ we are confident that the aforementioned results 
are simply due to the lowering of the deconfinement scale Tg as ^ is increased and not to the detailed analytical form 
of To(/i). Thus our picture should be qualitatively robust. 

Before closing this Section we make a brief comment on the possible study of the scenarios discussed above on 
lattice. Recently the density of states (DOS) method has been used to investigate the QCD phase transition at large 
fj, [13] ■ In this paper the QCD phase diagram is mapped by studying the plaquette expectation value in the /i — T 
plane. Although the lattice size implemented in [l3l | is relatively small and a finite volume study is still missing, so 
that the results should be taken as preliminary, an interesting phase transition is observed as fj, crosses a critical value 
fic at a fixed temperature. Moreover the quark number shows a sudden rise as fi reaches fic- The qualitative behavior 
is similar in the PNJL model, see Figs. [1] [5] and [31 In the PNJL calculation with parametrization of Case I and II 
(fixed values of Tq) at low temperature and in the correspondence of the the chiral crossover a small jump of the 
Polyakov loop occurs, the true Z{3) crossover being shifted to larger values of ^. On the other hand, in Case III with 
a ^-dependent Tq, a net rise of $ occurs in correspondence of the chiral crossover. It would be very interesting if by 
means of the DOS method one could compute the expectation value of the Polyakov loop, as well as the chiral and 
the Polyakov loop susceptibilities, in the low temperature regime as a function of /x. This lattice calculation might 
improve the understanding of the new low temperature/large chemical potential state of matter claimed in jl3s |, and 
at the same time it would allow to distinguish between the two PNJL scenarios discussed in this paper. 



VI. CONCLUSIONS 



In this paper we have investigated on the landscape of the possible phases of the neutral two flavor PNJL model. 
We have considered the logarithmic effective potential of the Polyakov loop U \2^, |25[ , see Eq. ^ , and a non local 
interaction in the quark sector, see Eqs. P^ - fH)) . Our main results are summarized in Figs. [3] and [T] Fig.|3]corresponds 
to a fixed value of Tq in the Polyakov loop effective potential. In this case the phase diagram is qualitatively similar 
to that obtained in the large Nc approximation of QCD [isj . 

In particular, at high chemical potential and low temperature we find a phase in which the main contribution to 
the thermal quark population is given by Z{5) neutral states, that is three quark states made of one red quark, one 
green quark and one blue quark. This characteristic resembles the quarkyonic phase of Ref. [45| . The quarkyonic-like 
structure of the ground state of the PNJL model has been already noticed in Ref. ^331] in non neutral and three flavor 
version of the model. Moreover the Z{3) transition line has been already studied in Ref. [s^l with a different effective 
potential for the Polyakov loop and in a non neutral state. The results of Ref. [SO^ are qualitatively similar to ours. 
Therefore we suggest that the quarkyonic-like state of matter is a feature of the PNJL model, independently to the 
number of flavors and to the difference of the chemical potentials between quarks, as far as a /i dependence of the 
coefficients of U is not considered. 

In Fig. [7] we show the phase diagram of the model when a /i dependence of the coefficients of the effective potential 
of the Polyakov loop is introduced. We have used the analytic form suggested in Ref. [1^ . The main results are the 
lowering of the Z{3) transition line of Fig. [H and the shrinking of the quarkyonic-like phase window of Fig. U) We 
have used the form of To(/i) of Ref. [11]. We believe that the result is rather robust as it does not follow from such a 
detailed form but only from the lower scale of deconfinement when mu increases. 

We have not considered in this work for simplicity the possibility of color superconductivity at high /i [ssl . [s^ . 
At a first sight it could seem that the results found with parameterizations I and II, i.e. a quarkyonic-like phase 
at high chemical potential and low temperature, exclude the possibility of a superconductive gap in the spectrum. 
This reasoning could be supported by the observation that the quarkyonic-like phase is similar to a confined phase, 
differing from the latter only for a non zero value of the quark density. Such a conclusion is not necessarily true. As a 
matter of fact, even if not noticed explicitly in Ref. [1^ for the two flavor and in [stI, [i^ for the three flavor models 
where the color superconductivity has been kept into account, in the quarkyonic-like region (high /i and small T) the 
minimization of the thermodynamic potential leads to a phase where quarks have a color superconductive gap in the 
spectrum. It is the 2SC gap [SSt in the two flavor case, and the CFL gap [s^ in the three flavor case. Therefore the 
realization of a color superconductive phase in the PNJL models at high jj, and small T is not forbidden in principle. 
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even if the ground state has a quarkyonic structure. 

An interesting investigation is the computation of the spectra of the mesonic and baryonic thermal excitations in 
the quarkyonic-hke phase of the PNJL model, and compare them with those obtained in a different model [i^ that 
mimics QCD in the large Nc approximation. We are now working on this topic and the results will be the object of a 
forthcoming paper. 
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